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ABSTRACT 


The  behaviour  of  a  surface  density  wave  situated  in 
parallel  electric  and  magnetic  fields  in  a  homogeneous  plasma 
of  a  cylindrical  intrinsic  semiconductor  bar  is  studied.  The 
possibility  of  finite  plasma  densities  at  the  surface  of  the 
semiconductor  is  included.  A  particularly  simple  model  is  used, 
requiring  only  a  thermal  equilibrium  electron-hole  plasma  and 
low  recombination  surface.  Explicit  solutions  of  the  basic 
equations  are  presented.  The  dynamic  behaviour  of  the  plasma 
is  investigated  by  perturbation  theory.  The  critical  field 
at  which  instability  appears  and  the  frequency  of  the  helix 
which  thereby  occurs  are  determined.  Helical  waves  arise  when 
the  field  exceeds  a  critical  value.  In  this  calculation,  linear 
terms  in  the  induced  azimuthal  magnetic  field  are  retained. 

Two  extreme  cases  of  reasonably  low  and  sufficiently  high 
applied  electric  field  are  considered. 
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NOTATION 


R  radius  of  semiconductor  sample, 

n  electron-hole  density. 

n  density  at  which  no  net  recombinations  or  generations 

take  place. 

n  density  at  the  boundary  at  which  surface  recombin- 

s 

ations  do  not  take  place. 

£  rate  of  creation  or  recombination  of  electron-hole  pairs, 

e  absolute  value  of  electronic  charge. 

T  particle  flux. 

E  electric  field. 

5  magnetic  field. 

A,$  parameters  associated  with  azimuthal  magnetic  field. 

T+  temperatures  of  holes  and  electrons. 

y+  mobilities  of  holes  and  electrons. 

D+  coefficient  of  diffusion  of  holes  and  electrons . 

s  surface  recombination  velocity. 

v  =  v+  -  v”  the  sum  of  the  absolute  values  of  the 

z  z  z 

longitudinal  velocity,  v  is  positive  for  positive 
electric  field  and  is  proportional  to  the  current. 
n&  average  electron-hole  density. 

Io  modified  zeroth  order  Bessel  function. 

Ji  first  order  Bessel  function. 

Bo, Bi  first  and  second  zeros  of  Io . 

Ko  first  zero  of  Ji . 
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n  tt 


velocity  of  light. 

velocity  of  sound  in  semiconductor. 

=  o)i  +  io)2  frequency  of  the  helical  wave, 
wave  number  of  helical  wave  in  azimuthal  direction, 
wave  number  of  helical  wave  in  z  direction, 
and  "+"  refer  to  electrons  and  holes. 
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Chapter  1:  Introduction 


1.1  Description  of  the  system 

The  system  to  be  studied  consists  of  an  electron- 
hole  plasma  situated  in  parallel  electric  and  magnetic  fields. 
The  plasma  is  contained  in  a  uniform  cylindrical  intrinsic 
semiconductor  with  diameter  much  smaller  than  the  length  of 
the  sample.  The  diameter  of  the  sample  is  of  the  order  of  a 
millimeter  and  the  length  is  of  the  order  of  a  centimeter,  so 
that  the  end  effects  may  be  neglected.  The  electric  field  is 
produced  by  electrodes  at  the  end  of  the  semiconductor  bar 
and  the  magnetic  field  by  an  electromagnet  which  is  rotatable 
in  a  horizontal  plane  to  complete  the  alignment  of  sample  and 
field  or  by  solenoid  of  diameter  barely  larger  than  that  of  the 
sample  with  which  used. 

When  a  charged  particle  moves  through  a  region  where 
an  electric  field  and  a  magnetic  field  are  present,  the  particle 
is  subject  to  two  forces.  The  electric  force  is  parallel 
to  the  field.  The  magnetic  force  is  at  right  angles  both  to  the 
velocity  of  the  particle  and  the  magnetic  field  strength.  When 
the  electric  field  and  magnetic  field  are  simultaneously  present 
and  parallel  to  each  other,  the  two  motions  are  combined  and  the 
final  path  will  be  a  helix.  Whether  the  helix  is  right-handed 
or  left-handed  depends  on  the  sign  of  the  charge  . 
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In  a  plasma,  the  longitudinal  electric  field  tends 
to  separate  the  superimposed  positive  and  negative  screws 
axially,  which  is  equivalent  to  a  rotation  of  one  screw  relative 
to  the  other.  The  screws  of  opposite  charges  will  attenuate  each 
other.  The  resulting  charge  separation  creates  both  an  azimuthal 
and  a  radial  electric  field.  In  the  linear  approximation,  these 
fields  together  with  the  longitudinal  field  act  on  the  unper¬ 
turbed  distribution  to  produce  a  radial  and  azimuthal  flow 
of  particles.  If  the  unperturbed  distribution  has  a  radial 
gradient,  this  flow  can  feed  particles  from  the  main  distribut¬ 
ion  into  the  screw,  thus  producing  a  growth  of  the  perturbat¬ 
ion.  When  the  influx  of  particles  is  sufficiently  large  to 
overcome  the  dissipative  effects  of  diffusion  and  recombination, 
the  helical  wave  will  grow. 

It  is  expected  that  for  a  given  value  of  electric 
field  there  will  be  a  threshold  value  of  magnetic  field  above 
which  growth  is  observed.  The  effect  corresponding  to  the 
steady  state  solution  is  given  in  Chapter  3-  As  the  applied 
longitudinal  magnetic  field  is  increased,  a  screw  type  state 
is  observed.  This  helical  state  will  be  unstable  for  suffic¬ 
iently  large  magnetic  field.  This  helical  state  solution  is 
given  in  Chapter  4. 

The  variables  of  this  problem  are  the  density  of 
electrons  and  holes,  the  electric  field  and  the  magnetic  field. 
The  electric  and  magnetic  fields  are  externally  applied  and 
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are  treated  as  independent  variables.  The  density,  the  internal 
electric  field  and  the  induced  magnetic  field  are  the  dependent 
variables . 

1.2  Historical  background 

The  helical  plasma  density  wave  (instabilities)  was  first 
proposed  by  Kadomtsev  and  Nedospasov1  as  an  explanation  of  a 
gas  discharge  in  a  longitudinal  magnetic  field  for  the  anomalous 
diffusion  in  the  positive  column.  Johnson  and  Jerde2  have  given 
this  theory  a  rigorous  mathematical  foundation.  Glicksman3 
suggested  that  a  slightly  generalized  form  of  Kadomstev  and 
Nedospasov  theory  could  explain  the  oscillist or ' effect  in 
semiconductors.  In  its  simplest  form  the  oscillistor  is  a 
semiconductor  bar  that  exhibits  thermal  voltage  and  current 
oscillations,  when  an  electron-hole  plasma  is  created  by 
injection  or  light  generation  or  collision  ionization  and 
sufficiently  large  electric  and  magnetic  fields  are  applied. 
Current  oscillation  has  been  observed  by  several  authors4*5 
when  a  semiconductor  plasma  is  placed  in  electric  and  magnetic 
fields.  The  helical  motion  of  the  plasma  has  been  experiment¬ 
ally  confirmed  by  Okamoto  et  al6  and  Misawa  et  al7.  Glicksman 
treated  only  the  case  of  an  injected  plasma  in  the  absence  of 
a  thermal  background.  However,  the  background  plasma  may  be 
of  great  importance  in  a  theory  of  the  helical  instability  in 
a  semiconductor  plasma,  as  has  been  pointed  out  by  others6*7. 
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Holter8  applied  Johnson’s  method  and  developed  a  theory  for 
helical  instabilities  in  solid  state  plasma  that  cover  all 
injection  levels  in  intrinsic,  n-type  and  p-type  semiconductors 
as  well  as  insulators  by  taking  the  background  of  thermal 
carriers  into  account.  A  quantitative  mathematical  description  of 
the  helical  wave  in  the  surface  density  mode  was  formulated 
by  Hurwitz  and  McWhorter9.  Recently,  the  case  of  an  inhomo¬ 
geneous  injected  plasma  in  a  cylindrical  Ge  rod  possessing 
low  surface  recombination  has  been  discussed  by  Schulz10. 

According  to  the  type  of  wave  investigated,  the 
existing  theories  can  be  separated  into  two  groups.  The  theory 
introduced  by  Glicksman  transferred  the  theory  originally 
developed  by  Kadomtsev  and  Nedospasov  for  a  helical  density 
wave  in  the  positive  column  of  a  gas  plasma  to  electron-hole 
plasma  in  a- semiconductor .  It  is  assumed  that  recombination 
is  dominant  at  the  surface  through  the  boundary  condition  that 
the  density  at  the  surface  is  zero  and  only  bulk  generation 
is  taken  into  account.  Modification  is  extended  by  other8, 
investigates  a  so-called  gradient  density  wave  which  needs  a 
gradient  in  the  steady  state  carrier  distribution,  usually 
caused  by  a  high  surface  recombination  of  the  crystal.  The  second 
theory  is  introduced  by  Hurwitz  and  McWhorter  for  a  homogeneous 
plasma  of  a  cylindrical  crystal.  It  is  assumed  that  electron-hole 
recombination  and  regeneration  depend  on  the  local  density. 

If  the  local  density  is  larger  (smaller)  than  some  value,  net 
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recombination  (regeneration)  will  take  place.  A  low  surface 
recombination  velocity  is  taken  into  consideration  through 
the  boundary  condition  at  the  walls. 

None  of  the  theoretical  works  in  this  field  has 
included  the  forces  due  to  the  internal  magnetic  field  in 
the  azimuthal  direction.  It  is  assumed  that  the  applied 
longitudinal  magnetic  field  is  usually  very  large  compared 
with  the  azimuthal  field  and  the  latter  may  be  neglected. 

As  a  consequence  of  the  omission  of  this  force  ,  we  can 
expect  that  the  results  of  the  above  authors  would  be  modified 
for  large  electric  current  when  the  induced  azimuthal  field 
becomes  comparable  to  the  longitudinal  field. 

1.3  Outline  of  the  method 

In  the  present  calculation,  we  treat  the  surface 
density  wave  by  including  the  induced  magnetic  field.  The 
induced  magnetic  field  is  assumed  to  be  small  and  only  linear 
terms  are  kept.  The  behaviour  of  the  plasma  is  investigated 
by  perturbation  and  the  effects  of  the  induced  magnetic  field  are 
found  in  the  calculation.  We  follow  Hurwitz  and  McWhorter  by 
considering  recombination  and  regeneration  and  assuming  a 
small  surface  recombination  velocity.  Hurwitz  and  McWhorter 
assumed  that  the  boundary  properties  are  not  different  from 
the  bulk.  This  approximation  is  inappropriate , for  the  density 
at  the  surface  may  be  different  from  that  in  the  bulk.  We  will 
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consider  that  there  is  a  difference  in  bulk  and  surface  den¬ 
sity.  Thus,  two  boundary  conditions  are  to  be  applied.  The 
first  consider  relates  the  flux  at  the  boundary  with  the 
surface  recombination  velocity  and  the  second  condition  ensures 
that  the  total  rate  of  generation  of  electron-hole  pairs  is 
equal  to  the  total  rate  of  recombination. 

The  mode  we  are  considering  involves  only  the  thermal 
equilibrium  electron-hole  plasma.  The  mobilities  of  electrons 
and  holes  are  considered  to  be  determined  by  collision  of  free 
electrons  and  holes  with  the  lattice  vibrations  alone  and 
collisions  between  electrons  and  holes  are  neglected.  When  a 
drift  velocity  is  much  less  than  the  velocity  of  sound  in  the 
medium,  the  temperatures  of  electrons  and  holes  are  nearly 
equal  to  the  temperature  of  the  lattice.  The  mobilities  and 
temperatures  of  electrons  and  holes  are  approximately  constant 
at  a  reasonably  small  electric  field.  As  the  electric  field  is 
increased,  the  electrons  and  holes  may  then  be  considered  as 
becoming  hot  and  out  of  thermal  equilibrium  with  the  lattice. 

A  significant  departure  of  mobilities  and  temperatures  from 
their  low  field  values  are  then  expected.  The  mobilities  will 
now  vary  inversely  with  the  square  root  of  the  electric  field 
and  the  temperatures  vary  proportionally  to  the  electric  field1 

To  solve  the  problem,  we  make  the  assumption  of  quasi 
neutrality.  The  differences  in  the  density  fluctuations  of  the 
two  types  of  charges  will  be  much  smaller  than  the  fluctuations 
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themselves.  We  may  neglect  the  existence  of  space  charge  with¬ 
out  excluding  the  internal  electric  field,  for  a  very  small 
space  charge  may  set  up  a  large  internal  electric  field.  These 
fields  are  established- in  a  direction  such  that  they  try  to 
equalise  the  densities.  In  the  present  problem,  we  want  to 
investigate  the  density  fluctuations  that  are  equal  for  both 
types  of  charges.  We  neglect  the  difference  in  the  density  of 
the  charges  (i.e.  no  space  charge)  and  still  assume  that  an 
internal  electric  field  exists.  This  assumption  is  known  as 
quasi-neutrality . 

The  equations  of  motion  of  the  system  are  obtained 
from  Euler’s  equation  of  hydrodynamics.  We  treat  the  plasma 
as  two  fluids  and  write  an  equation  for  each  of  the  two  fluids, 
electrons  and  holes.  The  other  hydrodynamics  equation  which 
we  must  use  is  the  continuity  equation  for  a  fluid.  Using  the 
equations  of  motion  and  the  equations  of  continuity  for  holes 
and  electrons,  the  dispersion  relation  governing  the  growth 
and  propagation  of  the  helical  wave  will  be  investigated. 

We  study  the  steady  state  and  the  time  dependent 
solutions  of  the  system.  In  the  steady  state,  the  density  at 
every  point  in  space  remains  constant  in  time.  We  assume  that 
density  fluctuation  can  be  written  in  the  form  of  a  series  of 
zeroth  order  Bessel  functions  and  we  only  retain  terms  linear 
in  self  induced  magnetic  field.  To  study  the  time  dependent 
solution,  we  assume  that  deviations  in  density  and  in  electric 
field  from  the  steady  state  values  are  small  and  we  only 
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retain  terms  linear  in  these  deviations.  We  assume  a  screw 
type  solution  for  the  time  dependent  part  of  the  densities 
and  use  the  perturbation  approach  to  find  the  effect  of  self 
induced  magnetic  field  on  the  densities.  Conditions  for  the 
instabilities  in  the  screw  type  motion  are  then  investigated. 
The  condition  for  stability  is  found  by  requiring  the  imaginary 
part  of  the  frequency  to  be  equal  to  or  greater  than  zero. 

This  equation  contains  the  wave  number  of  the  helical  wave, 
the  critical  longitudinal  magnetic  field,  the  critical  longit¬ 
udinal  electric  field  and  the  induced  azimuthal  magnetic  field 
as  parameters . 
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Chapter  2:  Basic  Equations 

We  consider  an  infinite  uniform  cylindrical  intrinsic 
semiconductor  (i.e.  total  number  of  free  electrons  and  holes 
are  equal)  of  radius  R  so  that  the  end  effects  may  be  neglected. 
We  assume  the  electrons  and  holes  behave  as  a  non-degenerate 
gas.  In  this  cylinder,  no  electrons  and  p0  holes  are  distri¬ 
buted  uniformly  in  thermal  equilibrium.  Since  the  semiconductor 
is  an  intrinsic  one,  we  may  write  n0  =  po  =  n  and  the  following 
condition  must  therefore  be  fulfilled, 

Vno  =  Vpo  =  Vn  =  0. 

When  an  electric  field  is  applied  to  the  semiconductor,  the 
thermal  equilibrium  is  upset.  Quite  large  changes  of  n  and  p 
may  however  be  obtained  by  applying  a  voltage  of  appropriate 
sign  to  the  electrodes.  A  quasi-neutral  density  perturbation 
is  created  under  the  influence  of  an  applied  electric  field. 
Electrons  and  holes  are  equal  in  number  only  in  intrinsic 
semiconductors,  while  after  excitation,  the  electron-hole  pairs 
produced  last  only  for  a  characteristic  time  called  the  life 
time.  Let  the  disturbed  electron  concentration  n  be  equal  to 
no  +  ni,  and  the  disturbed  hole  concentration  p  be  equal  to 
Po  +  Pi-  Since  we  assume  no  =  Po  =  n,  unless  ni  =  pi  a  space 
charge  p  will  be  set  up  giving  for  homogeneous  material 


p  =  e  ( p  i  -  n  i ) 
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We  assume  that  m  must  be  nearly  equal  to  pi  unless  a  very  strong 
electric  field  is  present.  This  can  only  occur  near  the  surface 
or  in  a  region  in  which  the  impurity  concentration  varies  very 
rapidly ,  giving  rise  to  a  high  concentration  gradient.  We  have 
from  Poisson’s  equation 


di v  E  =  j  =  47re(pi  -  m)  •»  , 

where  e  is  the  permittivity. 

For  the  present,  we  assume  the  generated  plasma  created 
by  the  applied  electric  field  is  quasi-neutral.  Under  the  assum¬ 
ption  that  the  times  and  distances  of  interest  are  much  larger 
than  the  mean  free  times  and  paths  of  the  free  carriers ,  we  con¬ 
sider  the  behaviour  of  a  semiconductor  in  which  a  deviation  from 
the  equilibrium  carrier  concentration  has.  been  created.  Clearly, 
if  the  influence  causing  the  deviation  is  removed,  the  concent¬ 
rations  will  return  after  a  time  to  their  equilibrium  values. 

Let  the  rate  of  creation  or  recombination  of  the 
electron-hole  pairs  be  E,.  We  then  expressed  the  rate  of  change 
of  n  and  p  by  mean  of  the  equation12 

d(n-n) 


dt 


=  -(n  -  n) E  j 


Re com . 


where  n  is  the  electron-hole  density, and  n  is  the  density  at 
which  no  net  recombinations  or  generations  take  place.  The  term 
-(n  -  n)E  takes  account  of  creation  and  recombination  of 
electron-hole  pairs.  We  assume  no  trapping  of  electrons  and 


holes  so  that  electrons  and  holes  are  generated  in  pairs. 
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The  continuity  equations  for  electrons  and  holes  are 


(1) 


The  subscript  and  "+"  denote  the  electrons  and  holes  res- 

pectively.  r+=  nv+  are  the  particle  flux  vectors  and  v+  is 
the  average  velocity  of  either  component  of  the  plasma. 


The  equations  of  motion  are 


?+  +  D+Vn  +  y+n8  +  y  +  f  +  x  8=0 


(2) 


where  we  have  used  the  quasi-neutral  assumption  and  assumed 
that  lattice  scattering  is  dominant  and  neglected  the  collis¬ 


ions  between  electrons  and  holes.  8  and  8  are  the  electric 


and  magnetic  fields.  y  +  are  the  mobilities  and  D+  are  the 
diffusion  coefficients  of  electrons  and  holes.  y+  and  D+  are 
related  by  D+  =  (KT  +  /'e)y+J  which  is  known  as  Einstein's  relat¬ 
ionship.  T+  are  the  temperatures  of  electrons  and  holes. 

Solve  equations  (2)  for  f+  explicitly  and  substitute 
in  equations  (1),  the  continuity  equations  are  then  written 


V  •  [  y +D+8( Vn *8)  +  y^y|n8(8»E)  +  DjVn 

+  D|y  +  (8xVn)  +  y_J.nET  +  y+y|n(BxE)  ] 

=  _  +  (n-n)5  , 


(3) 


where 


1  +  y2B2 


1  +  y 2B 2 


Vi 


and 
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We  have  used  cylindrical  coordinates  for  reasons  of  symmetry  in 

solving  ?+ .  The  z  axis  is  taken  along  the  direction  of  applied 

electric  and  magnetic  fields.  We  have  neglected  [Bfl(R)/Bo  ]2 

u  z 

in  comparsion  with  unity.  Here  Bg(R)  is  the  induced  azimuthal 
magnetic  field  at  the  boundary  R.  is  simply  related  to  the 
electric  current  through 


e  fR 

BQ(r)  =  -  n(r’)v  rT  drT  , 

H  2  n  2 

c  r  J  0 

^  -[•  ->■—  i 

where  v  =  v  -  v  and  v  are  the  z  components  of  the  drift 
z  z  z  z 

velocities  and  assumed  to  be  independent  of  r  and  z.  In  the 

literature  Bg(r)  has  usually  been  neglected.  We  assume  a  simple 

form  of  BQ(r) 
u 


where  n  is 
a 

BQ(r)  seems 
y 


e  v  n  r 

Bfl(r)  =  - 5 — - —  ,  (4) 

0  c2 

the  average  density  of  charges .  The  above  form 
reasonable  if  v is  assumed  independent  of  position. 


9  Irffi 
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Chapter  3:  Steady  State  Solution 


In  this  case,  there  are  no  variations  in  time  and 
in  space.  The  plasma  will  occupy  the  available  volume  with  a 
distribution  in  density  which  depends  on  the  sources  and  sinks 
for  the  plasma.  The  sink  of  the  plasma  will  be  primarily  at 
the  surface,  with  some  volume  recombination  also  allowable, 
so  that  the  plasma  density  will  have  a  maximum  at  the  center 
of  the  cylinder  and  fall  off  to  some  smaller  value  at  the 
surface.  We  assume  that  in  steady  state  Eg=0,-^-=0,^-=0 


and  ^ -  =  0.  Where  cylindrical  coordinates  are  chosen  with  the 


9  z 


axis  along  E0  . 

z 

In  the  steady  state  the  following  condition  must  be 

satisfied . 

(a)  We  neglect  the  induced  azimuthal  magnetic  field  and  write 
the  r  components  of  equations  (3),  we  obtain 


By  eliminating  Eor  between  the  two  equations  of  (5),  the  result¬ 
ing  equation  is 


where 


(6) 


Solution  of  (6)  is  of  the  form 


.  f.  0  5'lfi  89vt£f:  O'fQOD  I.soirrbfliJ  3rr®: 


. 

■ 


ax  no 


•3  t  r  • 
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I  1 )  r  ->  . 


14 


n  =  N0I o ( 3  or)  +  n. 


(7) 


where  N0  is  the  deviation  of  n  at  r  =  0 . 

To  complete  the  solution,  the  boundary  conditions 


?+(R)  =  s  [n (R)  -  n  ]  , 


(8) 


f2-n  f  R 

(n-n)  £rdrd0  +  2ttRs  [n ( R) -n  ]  =  0 

J  oh-  s 


(9) 


are  applied.  Where  r  +  (R)  are  respectively  the  radial  flux 

components  of  electrons  and  holes  at  the  surface,  s  is  the 

surface  recombination  velocity  and  n(R)-n  is  the  excess 

s 

carrier  density  at  the  surface. 

n  and  n  are  connected  by  the  first  boundary  con- 
s 

dition  through 


D  ’  -  D  ’  <a  -  \ 

-n„(d:  +  u;  „t  x  „T)  UJ-ofg-!Lr) 


±  y;  +  u 


dr 


r=R 


=  N0s10(3oR)  +  s(n-n  ) 

o 


(10) 


and  No  is  given  by  second  boundary  condition 


No  = 


-R|(g-ns) 


R 


0 


(11) 


Io(3or)rdr  +  R  I  <^(  3 aR) 


The  radial  component  of  the  internal  electric  field  E0r  niay 
be  determined  from  the  condition  that  in  the  steady  state. 


the  radial  flux  of  both  the  electron  and  hole  must  be  equal. 

E0. 


D+  D-  1  dn 


r  n  dr  • 


(12) 


- 

-■Da  +  (flof')olso  -  <T  -rr-  u  i«  +  ;  n  ■ .  Ti- 

■ 

no 
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(b)  Including 
components  of 


the  induced  azimuthal  magnetic  field,  the  r 
equations  (13)  are 


-(n-n)£  +  D|(-^.  +  -  ~)  +  y!  -  |-(rnEo  ) 
± 2  r  3r  ±  r  3r  r 

d  V 


+  y±y±  r  lF(rnEozB0)  =  °‘  (13) 

Eliminate  E  between  the  two  equations  of  (13),  the  equation 
is 


where 


2-  +  (-  +  Ar)~  -  ( 3  o 
~  2  r  3r 

3r 


-  2A)n  +  3on  =  0, 


(14) 


A  = 


y (  u_+  y+) 

y  |D_|_  +  y^Dj 


$Eo 


$  = 


B 

R" 


$ 


is  the  induced  magnetic  field  at  the  boundary  r  =  R.  Since 
n  depends  only  on  r,  we  may  replace  the  partial  derivatives 
by  the  total  derivatives .  Introducing 


N  =  n 


n 


5 


the  differential  equation  is  then  written  as 


—  +  (i  +  Ar)^  -  (Si  -  2A)N  =  0  .  (15) 

dr2  r  dr 

If  we  put  A  =  0,  which  corresponds  to  zero  induced  magnetic 
field,  equation  (15)  will  then  reduce  to  equation  (6). 

Now  consider  equation  (15)  using  method  of  perturb¬ 
ation,  the  next  approximation  can  be  found.  We  may  assume 


,  5  •  «(;  -  If  -  •!  <-:A  +  ~  ■+  — 

B  ■  " *a  . 

. 


vJ/.O.'  9  t  C 
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the  solution  of  (15)  as  a  series  of  zeroth  order  Bessel  funct¬ 
ion.  The  solution  is  written  as 

N0  Cl o (Bor)  +  a 1 1 o ( 3 i r) ] , 

where  3o  and  3i  are  determined  by  the  boundary  conditions. 
Substitute  this  expression  in  (15),  multiply  by  I0(3ir)rdr 

and  integrate  over  r  from  zero  to  R.  Neglecting  the  higher 
order  terms  in  $  and  using  the  orthogonal  property  of  Bessel 
functions,  we  obtain  an  expression  for  ai  whose  linear  term 
in  A  is 


A 

« 

r 2 -1  Io(3ir)dr 

0 

rR 

(3?  -  38) 

I  o  (  3  i  r )  r  dr 

• 

0 

Hence  approximation  to  the  solution  of  equation  (.12)  may  be 

written  as  6o 

n  =  N0 [I o ( 3  or )  +  ail0(3ir)]  +  -  n.  (17) 

6§  -  2A 

Using  the  value  of  3o,  together  with  the  boundary  condition 
(8),  3i  can  be  determined  from  the  following  relation 

-{(Dj.  +  ay  | )  No  |^[Io(3or)  +  a  1 1 0  (  3  i  r )  ]r=R 


y^y|(y_  +  y,) 

+  — 1 - —  $RE0_[No  Uo(3oR)  +  ai  I  o  (  3  i  R) } 

y :  +  y| 


3? 


35  -  2 A 


n  -  n  1} 


s  N0  [ I o  ( 3  o  R )  +  ai  I  o  (  3 1  R)  ]  +  s[- 


38 


35  -  2 A 


n  - 


n  ] 

s  , 


(18) 


■ 


. 


,  :  (i  £ )  noXoF.i-T  s  1C  ro  :a  £oa  c  -t  .10 ^Butfxoiqq*  on,  ' 


' 

>  ';  +  ’■  0-0  ‘  ±C  ' 


+  fdji  '  -•  -  X  'I  o-'I 


17 


where 


a  = 


Dj  " 

yj  +  v'_  * 


No  Is  determined  by  the  second  boundary  condition  (9)  and  is 
given  as  follows 


( 


2  p2  _ 

-S - l)n  5—  -  r|( 

2  o  rt  2  S  o  2 


3  0 


Nn  = 


3o  -  2A 


35  -  2A 


n  -n  ) 
s 


R 


•R 


Io(3or)rdr  +  ai  Io(3ir)rdr  +  R-p[I0(3oR)  +  aiIo(3iR)] 

-  O'  Jo  * 

(19) 

The  radial  component  of  the  internal  electric  field  is  obtained 
from  the  condition  of  steady  state. 

a  dn 


E0  =  —  -t—  -  ar$E0  , 
u  r  n  dr  z  3 


(20) 


where 


a  = 


y_y^  -  y+y{ 
vl  +  yj 


•JJ  + 

JD 
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Chapter  4:  Helical  Solution 

It  is  expected  that  for  a  given  electric  field  there 
should  be  a  critical  magnetic  field  and  a  critical  frequency 
at  which  growth  is  first  observed.  At  this  point,  which  we  shall 
call  threshold,  the  frequency  co  and  propagation  constant  k 
are  both  real,  corresponding  to  neither  growth  nor  attenuat¬ 
ion.  For  longitudinal  magnetic  field  below  the  critical  field, 
the  helical  mode  decays .  Above  the  critical  field  the  helical 
state  is  unstable.  The  imaginary  part  of  the  frequency  depends 
critically  on  the  longitudinal  magnetic  field.  In  this  calcu¬ 
lation,  we  shall  include  the  induced  azimuthal  magnetic  field 
up  to  the  lowest  order. 

We  begin  with  the  same  basic  equations;  the  equations 

of  motion 

f  +  +  D+Vn  +  y+nE  +  y  +  ?+x  ^  =  0  (21) 

and  the  equations  of  continuity 

+  y?±  =  -(n-n)5.  (22) 

From  equations  (21),  f+  are  obtained  explicitly,  and 
?+  may  be  substituted  into  the  continuity  equations.  To  obtain 
r+ ,  we  write  equations  (21)  explicitly  as  three  equations  in 
the  cylindrical  components  of  T+,  and  solve  for  the  components 
of  ?  respectively.  From  the  expressions  for  these  components, 
we  write  a  vectorial  expression  for  f+. 


This  expression  may 
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be  written  as 

r±  =  -  [  y2D|£(Vn«B)  +  y2y;nS(E^E)  +  DjVn 

+  D|y  +  (§xVn)  +  y^.nE  +  y+yjn(BxE)  ]  (23) 
Substituting  expression  (23)  in  equations  (22),  we  obtain 
V  *  C  y2D{S(Vn*§)  +  y2y|n§(§»S)  +  DjVn 

+  D|y+(i^xVn)  +  yjnS  +  yhy|n(5xS)  ] 

=  §£  +  (n-n) £  .  (24) 

Since  we  assume  the  induced  azimuthal  field  is  normally  much 
smaller  than  the  applied  longitudinal  field,  we  may  replace 
1  +  y2B2  by  1  +  y2B2z,  and  write 

M±  D± 

yt  =  -  an(}  D  '  =  - 

1  +  t+Bo2z  -  1  +  y2Bo2z 

We  obtain  from  (24) 

-  |^r  +  (n-n)£  +  D|V2n  +  y|V*(nE)  +  y+y|V*(nBxE) 

+  y  *D_J_B  •  V ( B  •  Vn )  +  y2y|B  •  V  (nB  *E)  +  y  +  D|V*(BxVn) 

=  0.  (25) 

To  treat  this  problem,  we  suppose  a  fluctuation  in  the  distri¬ 
bution  of  carriers  to  have  originated  at  a  certain  moment  of 
time,  and  introduce  the  perturbation  terms  ni  and  Ex  defined 

by 

n  =  n0  +  ni , 


E  -  E  o  +  5 1  . 


(26) 


* 


' 
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Where  nQ  and  20  are  the  steady  state  values  of  n  and  S,  and 
ni  and  Ei  are  their  time  dependent  parts,  and  are  assumed 
small.  Replacing  n  and  E  with  the  above  expressions  in  (25), 
and  from  which,  subtracting  the  steady  state  equation,  we 
obtain  to  first  order  of  perturbations  of  nj  and 

-ni€  -  +  D  ’  V2m  +  y'V-(n0£i+  niS0) 

+  y±y  +  V • [§x(n02i  +  ni20)]  +  y+D^§* V(S«Vni ) 


+  y +  y *  V  [B  •  ( n  0E  i  +  n  iS0 )  ]  +  y±D;v(SxVm)  =  0.  (27) 


To  obtain  the  normal  mode  for  the  problem,  we  now 
assume  a  helical  form  for  the  density  ni  and  the  electrostatic 
potential  Vi  from  which  the  field  Ei  is  derived.  The  general 
forms  in  cylindrical  coordinates  (r,6,z)  are 

ni  =  [C  +  f(r)]  exp[i(oot  +  m9  +  kz)] 

(28) 

and  Vi  =  g(r)  exp[i(tot  +  m0  +  kz)], 

where  a)  is  the  frequency  of  the  perturbation,  m  is  the  wave 
number  in  the  azimuthal  direction,  k  is  the  wave  number  along 
the  z  axis,  and  C  can  be  determined  by  the  boundary  condition 
( 9 )  at  the  walls  . 


-2 


C  = 


R 


0 


f(r)rdr  -  2R|  f(R) 


(29) 


R2  +  2R?- 


The  radial  part  of  ni  takes  the  form  f(r)  =  fo Jm(Kor) ,  where 
Jm(K0r)  is  the  Bessel  function  of  the  first  kind  of 
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order  m,  and  f0  is  an  arbitrary  constant.  The  electric  field  E 
may  be  expressed  in  terms  of  the  potential  V  according  to  the 
following  set  of  equations 

■it?  . 

im  _g  x 
r  r  » 


-  ik  |£  ex 
9r  > 


Ei  = 

A  ->o 

1 

.  iS  ex 

9r  5 

,  3E  i 

1  r 

r  9  0 

3E  i 

r 

92g  x 

3E  i 

r 

9r 

,  c  » 

9r 2 

3  z 

Ele  " 

.  Vi 

-  im  7  , 

i  3Ele 

r  90 

3(rEi 

9r 

0  ^  .  9g  x 

—  =  -im  r*  e  . 

3r  5 

3Ele  _ 

9  z 

E>z  ■ 

-  ik  Vi  , 

3Elz 

30 

3E,z 

9r 

.  i  9g  x 
=  -  ik  e  , 

3r  s 

9E  i 

z 

3  z 

X 

e  = 

exp[i(cat  +  m0  +  kz)]. 

2  g  X 

m  —  e 


x 


(30) 


Replace  BQ  by  its  lowest  order  expression,  that  is 

u 


the  term  linear  in  r 


B 


n  v  er 
a  z 


and  define 


B, 


$  = 


n  v  e 
a  z 


and  $ , 


We  obtain  from  (27)  by  keeping  only  linear  terms  in  B0 


z 


? 
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1  d  f  d(C+f(r)) 
r  dr  dr 


]  +  y; 


D±k2  +  iy ±kE o  z 


+ 


y±y jrE 


d( C+f ( r ) ) 
z  dr 


$  + 


j-  |?[r(C+f(r))Eor]  +  C  -5  -iu 

B.z 

+  iy±u>  -jX.  E«r  ][C+f (r)  ] 

[  +iy|y +m  BqzEoz 


+  2y^(y±Eoz±ikD±)  -  iky^rEo^  ][C+f(r)]  $ 


=  +  y±  f  aF(rn°af)  +  c  ±y±no 


m- 


B 


±  y+n0k  -  iy  +  y_J_m 


_z  drip 
1  dr 


]g(r) 


+  y +y j  [  2in0k  +  ikr  |~]g(r)$  ,  (31) 

where  E0r  has  been  derived  previously  in  equation  (20).  Since 
all  variables  depend  only  on  r,  we  have  replaced  the  partial 
derivatives  by  total  derivatives. 

We  substitute  Eor  in  equations  (31)  and  define  a  new 
variable  Z(r) 3  such  that 


g(r)  = 


1 

n0  (r ) 


&(r)  . 


(32) 


Rewriting  equations  (31)  in  terms  of  f(r)  and  £(r)  ,  we  obtain 

D !  —  [r-(-Cf£(-)^-]  +  u!  -  [a  -  ^-(C+f(r))] 

±  r  dr  ±  r  dr  n0dr 


+  [  -£  -  ioo  -  D|  —  -  D  k z  +  iky.Eo 


2  T 


+ 


Bo  -i  j 

+  iay.y'm  — —  —  ][C+f(r)] 

r  n0dr 


.  I?  /  -i  ,  a  \  d  ( C+f  (r ) )  * 

+  y.ylrEo  (1  ±  £— )  — > — -= — - — —  $ 
z  y ,  dr 


‘  + 


+  [  +iy2y|mB0  E0  (1  ±  ^-) 
xx  x-  y  + 


+ 


2y±y|E 


(l 


' 
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±  2iky^D±  -  ikau±uj  ^  ][C+f(r)]  t 


=  + 


K 


i  d  (  an, 

r  dr  rdr 


1  d_/  r_ 
r  dr^n0 


dnQ 

dr 


)  Hr) 


+ 


1  drip 
no  dr 


]  &(r) 


+  y±y;[  2ik  +  ik  ]  £(r)  t  .  (33) 

Since  g(R)  which  is  essentially  the  amplitude  of  the  potential 
of  the  electric  field  at  the  boundary,  must  be  finite  and 
since  n(R)  has  finite  values,  due  to  the  boundary  condition 
at  the  surface,  £(R)  must  then  be  finite. 

We  may  write  the  solution  of  (33)  as  Bessel  function 
of  order  one 


f ( r )  =  f0  Ji(Kpr) 


(3^) 


A(r)  =  £p  Ji (K0r ) , 

where  fp  and  £p  are  independent  of  r .  We  substitute  (34)  in 
equations  (33),  multiply  by  Ji(K0r)rdr  and  integrate  over 
r  from  zero  to  R. 

Define 


P  = 


_  1  dnp  _  p 

-  - -  —  r  i 

no  dr 


+  P2AR‘ 


(35) 


where 


Pi  =  -  —  Rfr(n-n  ) 
n0  £  s 


dip ( 3  or ) 

dr 


1 

D  5 


P2 


_2 _ 

3p2R2 


R2  ,  „s  >  dip ( 3  o  r )  1 

~  +  ^  ““ d? -  D 


dip  (  gpr)  1  -I 
dr  D  J 


+  R^(n-ng)  a2 


3 


- 


'  J - ; 

o  s 


24 


and 

o  _  ai 

32  A  ’ 

d  =  [  i  +  ^  e§R2(|  +  Ij)  ][  |i  +  r|  ] 

We  obtain  two  algebraic  equations  whose  coefficients  are 
functions  ofthe  parameters  decribing  the  plasma  and  the  wave 
parameters  go  and  k  only. 

[  D|{M3  -  m2(M4  -  Ms)}  +  ay j { (M6  -  M8 )  +  (M7  -  M9)dR2$} 

+  C  -£  -ia)  -  D+k 2  +  ikp  Efl  )(Mi  -  M2)R2 

—  —  z 

+  iay+y|B0  m(M1Q  -  Mu)  ]  f0  +  C  {+imy2y;B0  E0  (1  ±  ^— ) 

+  2y  p|( 1  ±  ^ — )Eq  +  2iky'D+}(Mi  -  M2)R2 
-  -  y  +  ^  -  - 

+  U+U’E.U  ±  r->M:2R2  -  iaky.vi!(M,j  -  Mi„)R2  ]  f0  $ 

-  “  ^  “ 


=  [  +yj(M3  -  m2M4 )  ±  y;(M6  +  M7dR2$  )  +  y+k2MxR2 

-  imy +y|B0 zMi q  ]  £0  +  y  +  y|  [  2ikMiR2  +  ikMi 3R2  ]  Z0  $  . 

136) 

In  equation  (36) 


y  +  y 


d  = 


Eo 


— (T  +  T  )  ’Z 

e  +  - 


and  all  the  MTs  are  expressed  as  integrals  of  Bessel  function. 


fR 

MiR2  = 

■  0 

fR 

J?(Kor)rdr  , 

m2r2  = 

lo 

CiJi(Kor)rdr  , 

- 


tK 

. 


'  • 


]  ^  o'*  [  (.  “  .  C;':iM+UX  t 


: 


. 


*  '"I  -  ]  = 


+  4 


- 

L  ' 


n  i  ..s  )h  '} 


«  TOlC  C>IS, 
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where 


M- 


fR 


0 


jb[rdJl^°r)]  J 1  (Kor ) dr  , 


Mi 


rR 


,o 


-  J i ( K0r ) dr  , 


M. 


rR  c 
f  £l 
-  0  r 


J i ( K0  r ) dr  a 


M, 


R 


i  0 


I^CrPiJ^Kor)]  J  i  (K0r  )dr  , 


M- 


R 


J  0 


d_ 

dr 


[rP2Ji(K0r)]  Ji(K0r)dr  , 


M, 


-R 


0 


d_ 

dr 


[ rC i P ! ]  J i (K0r )dr  , 


*R  d 


0 


dr 


[rCiP2]  Ji(Kor)dr  , 


Mi 


i  o 


rR 


0 


Pi  J i (K0r) dr  , 


M, 


rR 


0 


CiPi  J i ( K o r ) dr  , 


Mi  2R 


2  _ 


R  2  dJ  t (Kor )  T  fv  \  n 
r  — v  ■"  -  -  J  i  (K0r)dr  , 

0 


rR 

Mi 3R2  =  I  r2Pi  J?(K0r)dr  , 

J  0 


m14r2  = 


R 


0 


r2CiPi  Ji ( K0r ) dr  , 


C 


C  i  f  o  • 


(37) 


.  '  ■  .  ' 

. 

' 

.  •  : 

v  •  >  ' 

1 


. 


•  . 
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Dividing  equations  (36)  by  y_[_(M!-M2)  and  defining 


N  i 

M3  -  m2(M4  -  M5) 

Mi  -  M2  j 

n2 

M6  -  M8 

Mi  -  M2  j 

n3 

_  M10  -  Mu 

Mi  -  M2  5 

n4 

M  1  2 

Mi  -  M2  a 

N5 

_  M  i  3  -  M  i  4 

Mi  -  M2  ’ 

n6 

_  M  2  /  o  3 

Mx  -  M2  s 

n7 

M7  -  m9 

Mi  -  M2  j 

n8 

M8 

Mi  -  M2  5 

n9 

m9 

Mi  -  M2  J 

N  i  o 

m2M5 

Mi  -  M2  s 

Nn 

_  Mu 

Mi  -  M2  5 

N  1  2 

Mi  3 

Mi  -  M2  5 

we  obtain  two  equations  in  fo  and  £0  of  the  form 

A+  fo  +  B+  £o  =  0  } 


A  f  o  +  B  £  o 


0 


3 


(39) 


p9  cw~  s  ey 


■  if  w 
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where 

KT  r  .  y4  KT 

A  =  — ^  N i  +  aN2 - ^  R2  -  ^  R2  -  -*r  — 11  k2R2 

±  e  y;  y;  y.j.  e 

y  + 

+  r-r  ikE0_R2  +  imay  Bo  N3  +  [  +iy  2B  o  _E0  mR ( 1  ±  f—) 

M_j_  ^  z  —  z  z  y  ^ 

KT+ 

+  2)1  Eo  R(-l  ±  2_)  +  °ikRvi  — -  +  y  +  E0  R(1  ±  £-)N„ 

-  "  y+  xe  -z  y  + 

-  iakRy+N5  +  aN7Rd  ]  R  $,  (40) 

y  + 

B+  =  ± (N !  +  N i o )  +  (N2  +  Ns)  +  rn-  k2R2 ( 1  +  N6) 

^  ± 

+  imy +Bo  z ( N  3  +  N11)  +  [  -2iy+kR(l  +  Ne) 

-  iy+kR(Ns  +  N12)  +  (N7  +  N9)dR  ]  R  $. 

The  dispersion  relation  for  co  is  calculated  by  setting 
the  determinant  of  the  equations  equal  to  zero.  The  stability  of 
a  wave  is  investigated  by  treating  either  wave  number  or  the 
frequency  as  a  complex  quantity.  In  the  present  calculation, 
we  assume  that  the  wave  number  k  is  real  but  the  frequency  is 
complex,  thus  00  =  ui  1  +  iaj2. 

To  have  a  non-trivial  solution  for  f 0  and  &0 ,  we  set 
the  determinant  of  the  coefficients  of  equations  (39)  equal  to 
zero . 


(41) 


At  threshold,  tne  real  and  imaginary  parts  may  be 
separated.  We  then  obtain  two  equations  of  the  form 


•i^;r  ••  ]  +  i  ill  +  :  l)  o  ±4f-  + 

' 

•  J  •  -  5  >'  > 

- 


Oj  r)  '  i  >  *’  ' 
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A  i  o)i  +  B 1 0)2  +  C  i  i  =  0  j 
A20)i  +  B20)2  +  C22  =  0  s 


(42) 


where 


y  y 


A 1  = 


+ 


B2  =  (-r  -  — r)  C  m(N3  -  Nh)Bq  -  2kR(  1  +  N6  )  R  $ 


-  kR(Ns  +  N12)  R  $>]  , 


A  2  -  -  B 1  - 


-  [  -(Ni  +  NioMt^t  +  7Tt)  +  (N2  +  NsHttt  +  77-5-) 


y 


y :  y 


y  +  y. 


-k  2  R  2  ( 1  +  N6) 


2  tD  2 

1  .  1 


+  dR(N7  +  N9)(^y  +  7^-)  R  $  ]  , 


y 


KT,  '  KT, 

C11  «  [  —  N,  -aN 2  -  -S-R2-  K  — ik2R2  +  2u  E„  R(1  +  2_)  R  $ 

“  M_|_  “  T  &  M  _|_ 


+  y+E»zR(1  +  ^  N^R  4  -  4  ]  C  -(Nl  +  Nlo) 

y_ 

+  (N2  +  Ns)  +  —r  k2R2(l  +  Ns)  +  d(N7  +  N9)  R2  $] 

y  KT 

-  [  -  —  k  E0  R2  +  amy,B0  N3  -  2y  kR  -  R  $ 

y+  z  <  z  t  e 

-  ay  k  R2N5  $]  [  my  B0_(N3  +  N11) 

r  ■"  Zt 

-  2y_k  R2(l  +  N6)  $  -  y_kR2(N5  +  N12)  0] 

KT  p  y_ 

+  [  -  N 1  +  aN 2  -  -r  R2 - r  k2R2  +  2y  E0  (1  -  —  )R  $ 

e  yl  yl  -  z^  y_ 

+  y  E0  R2(l  -  — )  N4  $  +  aN7dR2  $  ]  [  -(Ni  +  Ni0) 

—  z  y 


i  k2R2 ( 1  +N6)  +  dR2 ( N 7  +  N9)  $] 


+  (N2+Ns)+ 


-■ 


' 


[  :  (eW  +  tVI)sH£  +  (aV'  +  I)  ^ 
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C22  = 


+ 


y_ 

[  — r  k  E0  R2  +  may  B0  N3 
y  '  z  -  u  z 


KT_ 

+  2y  kR2 -  $ 

e 


-  ay_kR2N5  $  ]  [  my+B0z(N3  +  Nil) 

-  2y+kR 2 ( 1  +  N6)  $  -  y+kR2(N5  +  Ni2)  <f>  ]  , 

KT,  r  y  ,  KT, 

[  — — N i  -  aN2 - j-  R2 - r  — — k 2R2  +  2y,E0  R2(l  +  — ) 

e  y|  yj  e  uz  y  ' 

+  y  E0  R(1  +  — )N4R  $  -  aN7dR^  0  ]  [  my  B0  (N3  +  Nn) 

'  z  y  —  z 

-  2y_kR2 ( 1  +  N6)  $  -  y_kR2(N5  +  Ni2)$  ] 

y,  KT 

+  [  -  —  kEozR2  +  may+BozN3  -  2y_|_kR2-^-  $ 

-  ay+kR2N5  $  ]  [  -  (Ni  -  N10)  +  (N2  +  N8) 

y 

+  —  k2R2 (  1  +  N6)  +  dR2(N7  +  N9)  $  ] 

M 

KT  p  KT 

-  [  — — N  i  +  aN  2  -  4-R2  -  4  - k2R2  +  2y  E0  R2(l-— ) 

e  y^  e  -  z  y_ 

+  y  E0  R2(  1  -  — )  N4  $  -  aN7dR2 $  ] 

-  uz  y- 

[  my  B  o  (N3  +  Nn)  -  2y,kR2(l  +  N6)  -  y ,  kR2  (Ns+Ni  2 )  O  ] 

+  z  +  + 

y_  2KT_ 

+  [  — r  kE0  R2  +  may  B0  N3  +  2y  kR  — —  $ 
y  z  ~  z  “S 

-  ay_kR2N5  $  ]  [  -(Ni  +  Nio)  +  (N2  +  N8) 
y, 

+  -4  k2R2 (  1  +  N6)  +  dR 2 ( N 7  +  N9)  $  ] 


(43) 


•  .  :  '  vi,  ,  *  !"  1  '-v.-  '' 


$  s  1 

■ 

' 

. 
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In  order  to  effect  a  considerable  simplification  in  the  mathe¬ 
matics,  it  is  assumed  from  here  onward  that  uiB2A  <  u2Bo 
'  ±0  ±  z 

<<  1,  we  may  replace  1  +  y2Boz  by  one,  and  all  the  primes  in 
(43)  may  be  dropped. 

From  equations  (42), the  imaginary  frequency  may  be 


written 


A 1 

1 

0 

*—• 

►— * 

CM 

< 

CM 

CM 

O 

1 

A 1 

B 1 

a2 

b2 

(44) 


The  threshold  for  the  onset  of  instability  will  occur  when 
(jo2  =  0  ,  because  the  to  variation  of  the  wave  first  changes 
from  attenuation  to  growth  at  this  point.  Therefore,  we 
set 


A i  C22  —  A2  C 


1  1 


=  0. 


(45) 


The  parameters  of  the  system  are  of  course  the  magnetic  and 

electric  fields  applied  to  the  plasma  and  the  induced 

azimuthal  magnetic  field.  Note  that  B0  now  represents 

z 

the  critical  magnetic  field. 

We  wish  now  to  calculate  both  the  extreme  cases 
for  weak  and  sufficiently  strong  electric  fields.  It  is 
known  that  the  electrons  *  and  holes ’ temperatures  and  mobilities 
are  functions  of  the  applied  electric  field  and  expressed 
for  very  weak  applied  electric  field,  as 


1 1 


V. 


Eo 


<<  c 


■ 


Jt  g  *3fiL  1rT  9fW 

'  ■  ■ 


S 

. 

sm<  f!?x9  e>rtd  r  xd  ;Ii/o  £  »  -  on  ri.  xw  9W 

>  8°  >>  =  c3 
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(46) 


for  sufficiently  strong  applied  electric  field  ,  as 


K  Eoz  » 


c 


s 


+  x  c  x  C‘  0^ 


(47) 


s  s 


Where  T°  and  y°  are  temperatures  and  mobilities  of  electrons 

and  holes  at  vanishing  electric  field,  c  is  the  velocity  of 

s 

sound  in  the  medium  and  T+  >  T°  .  In  weak  electric  field 

T  -  T°  is  proportional  to  Eo  »  thus  we  see  that  the  effect- 
—  —  z 

ive  temperatures  of  electrons  and  holes  are  raised  slightly 
above  the  equilibrium  value,  but  that  this  is  a  second  order 
effect,  so  that  the  usual  assumption  that  T+  =  T0  is  just¬ 
ified.  We  see  that  y+-  y°  are  also  proportional  to  Eoz  and 
that  provided  y°Eoz<<  c  ,  it  is  a  good  approximation  that  y+ 
is  independent  of  Eoz  and  equal  to  y°. 


Substituting  the  values  of  Ai,  kz  3  Ci i ,  C22  from 


(43)  in  (45)  and  after  performing  some  algebric  manipulation, 
we  obtain  for  both  weak  and  strong  applied  electric  field  the 
equation  for  the  requirement  of  stability. 


•  „  >r° 

' 

■ 

)!  ',  0  9  £=9  nT  .  J?  <  +T  brt£  j  n  sn  •  fi 

. 

s9lori  bn£  enoicfos 

- 

. 

. 

■ 

. 

si'"-'  <  iO  t$A  ei A  lo  ■I'L’.lsv  srlcf  gn  JuJliad  Z 
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b(y_+y+)BQz[  Hig-p  +  mH2Z^  +  H3Zlg^ 


+  H5Z?g|^  ]  +  H6  +  H7Z1  +  H 8 Z 1  +  H 9 Z 1 


3  /  2 

+  mH  4  Z 1  7  2 
=  0 


where 

7  _  k2R2 

Zl  ■  in 

Eo  R 

b  =  -"K"  - - 

N?  |(T+  +  T_) 

Hi  =  G?G4/  N1?  -  2G7GiGN^  , 

H2  =  -  GaGi  , 

H3  =  ~  [Gi(2G4G6  -  G 5 ) /N^  -  2G 7 ( G 1  +  G$G)Ni2] 
H4  =  G3G6  , 

H5  =  Ge (G4G6  —  G 5 ) /N 1  —  2G 7G 6N^  3 
H  6  =  G  2  G  j 

H7  =  -  Gi(2GeG  -  Gi)  , 

H 8  =  G  6 ( G  6  G  +  2G 1 )  j 
H9  =  -  Gi  , 

and 


G  =  1  -  3oR2/Ni  , 

Gi  =  ( N !  +  N10  -  N2  -  N8)/N!  , 
G2  =  N2/N!  , 

G  3  =  (N3  +  Nm)/Ni  , 

G4  =  2  +  N4  3 

G5  =  2(1  +  N6)  +  (N5  +  N12)  , 
Gg  =  1  +  N6  , 

G7  =  (N7  +  N  9 ) /N 1  j 


(48) 


(49) 


H<: 


, 


—  T 


fttaDrOS  -  f Vi\  (  e G 
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G  8  —  N  3 /N i 


if  y°Eoz  <<  C 


s 


b(y  +  +  ) B o  z 


EozR 


(y°  +  y ° )Bqz 


if  y  o  z  >>  c 


s 


b(y+  +  y_)Boz 


Equation  (48)  gives  the  smallest  value  of  E0  for  which 

z 

instability  may  occur.  The  general  behaviour  is  as  expected; 
larger  magnetic  fields  are  required  to  cause  instabilities 
when  lower  electric  fields  are  applied. 


At  the  onset  of  instability,  we  also  have  the  condit¬ 


ion  that  the  derivative  of  to2  =  0  with  respect  to  wave  number  k 
must  be  equal  to  zero.  When  this  condition  is  satisfied,  a 
large  number  of  waves  have  nearly  the  same  imaginary  parts, 
their  amplitudes  decay  or  increase  in  time  at  the  same  rate. 
Thus  producing  an  enhanced  effect .  Thus 


(50) 


Equation  (50)  together  with  equation  (48)  give  the  value  of 

k  =  k  in  terms  of  the  parameters  E0  and  Bo  ;  such  that  u>2 
c  z  z 

is  nearly  zero  for  the  neighbouring  points  of  k  . 


;  r  ■:  . 

. 
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Eliminating  E0  between  (48)  and  (50),  we  obtain  an 

z 

equation  relating  the  wave  number  in  z  direction  k  and  the 

critical  magnetic  field  Bo  and  the  induced  magnetic  field. 

z 

2H  5H  9  ^  Z !  9  7  2  +  31T1H4H9Z?  +  4H3H9  Zl7/2 

Boz  Boz 

+  m(H4H8  +  5H2H9 )  Z 1  +  ( 2H  3H  a  +  6H1H9  -  2H5H7)  ^  Zi^2 

Boz 

+  m(3H2H8  -  H4H7)Z?  +  4  ( H 1 H8  -  H5H6)  5^  z/72 

Boz 

+  m(H2H7  -  3H  4  H  6  )  Z 1  +  2(HiH7  -  H3H6)  ^  z/72 

Boz 

-  m  H2H6  =  0.  (51) 

Another  quantity  we  wish  to  calculate  as  a  function 
of  the  critical  magnetic  field  is  the  frequency  of  rotation  of 
the  helix.  We  solve  for  a)  1  from  equations  (42),  which  may  be 
written 


-C 1 1  B 1 

—  C  2  2  B  2 

mi  =  -  (52) 

A 1  B  1 

A  2  B  2 


which  gives  the  real  part  of  u)  at  the  stability  boundary.  Where 
Ai ,  A 2 ,  Bi,  B 2 ,  C 1 1  and  C22  are  given  in  (43)  and  noted  that 

Bi  =  -A  2  3 

B  2  =  A 1  . 

Therefore  a)  1  may  be  written 

A1C11  +  A2C22 

0)  1 R2  =  - 


A 2  +  A  2 


(53) 


' 


C  ~ 
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Substituting  the  values  from  (43),  after  much  algebra,  we 

obtain  an  equation  relating  the  real  frequency  uii,  the  electric 

field  (E0  )  variable  b  and  the  wave  number  (k)  variable  1\. 
z 

Here  again  we  consider  the  two  extreme  cases  for  weak  and  suff¬ 
iciently  strong  electric  field.  The  frequency  ooi  is  given  by 
if  u"Eoz  «  Cs, 


bi  [  Hi  i  +  Hi2b(y_  +  y+)<I>R  +  H13Zi  +  Hii+Zf  ] 

“  y_  “  y+)Boz  (  H i 5  +  Hi6Zi  +  H17Z1  ) 

+  |(u_  -  U+)  tR  (  Hi e  Z 1 1/2  +  H19Z,3/*  +  H2  0  2 1 5/2  ), 

(54) 


if  y+Eoz  >>  cs, 

bi  C  H11  +  Hi2b(y_  +  y+)$R  +  H13Z1  +  H14Z2  ] 


-  m(y  -  y  ) ( 1  +  — 
+  2 


y y+ 


y+  +  y 


— )B 0  (  His  +  H 1 6 Z 1  +  Hi7Zf  ) 

2  z 


+(y  -  y,)(i  + 

-  -t-  2 


y  y 


-K+  W-  ,  TT  „  V2 


■)$R  (  H 1 8  Z 1 


3/?  5 

/2  +  H20Z1 


+  H i  9  Z 1  '  2  "  u  17  /z 


), 


y+  +  y. 


(55) 


where 


b  1  = 


y  +  y 


y_y+ 


a)  iR‘ 


N 1  |(T+  +  T__) 


» 


H 1 1  =  Of  , 

H 1  2  =  — 2GiG7Nf  , 

H 1 3  =  "2G 1 G6  , 

H 1  4  =  G  6  ,  (56) 

Hi  5  =  G 1  [  Gs(l  -  G2)  “  GiGe]  , 


'  "  f'  ■;  f  ■ •  1 


'  1 16  ■  ■  '  6  %  T  .  i  V  '  #  1 10*3  3  8  '  if  ’ 

[  £:.  tSe  £.;  +  i>(  M  '+  a  dstH  +  nH]  id 

;  /H  +  i  S  a  x  H  +  e  1 1 '  "  _U  '\ 


- 


. 


r-b.:'  :  ( - 


4  1  -;M 


+  I ) ( 


+ 
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H  i  6  =  —  G3  (  G 1  +  G  6  —  G2G6  )  +  2GiGeG8  , 

H 1  7  =  G 6  (  G3  -  GeGe  )  , 

Hi s  =  Gi  [  Gi(2  +  Gg)  -  G 5 ( 1  -  G2)  ]  /  N f  , 

H 1 9  =  [  G 5 ( G 1  +  G6  -  G2G6)  -  2GiG6(2  +  Gg)  ]  /  , 

H20  =  G 6  [  G  6  ( 2  +  Gg)  -  G5  ]  /  Nf  . 


If  y°E0  <<  c  , 
±  z  s  * 


V*_  -  -  Vi  +  , 


y°  +  y° 


1  - 


y 0  y 0 


o)iR' 


-M+  2N 1  |  To 


if  y±E°z  >>  CsJ 


y-y+ 

(y  -  y  ) (1  +  - 

+  2 


y+  +  Vi. 


•)  =  0^  _  ^)d  + 


O  0^ 

y  y+ 


.  1  0  4.  0  TT1^ 

P+  +  P_  Eoz 


b  1  = 


_  Jk  (  32  A 


C  “  (  — — ) 

s  k  3ir 


p^  +  pf8 


u)  1 R  ‘ 


y0^y0^(y0  ^  ^  ^  ' 


Thus,  the  frequency  can  be  calculated  when  k  and  E0  are 

z 


determined . 
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Chapter  5:  Discussion 

Steady  state  solution  (17)  is  of  the  standard  form 

containing  Bessel  functions  in  the  limit  $  =  0.  However  the 

coefficient  N0  of  I0  vanishes  if  n  =  n  and  $  =  0.  If 

s 

$  =  0 ,  the  steady  state  solution  gives  an  increase  in  density 
on  the  axis  and  a  decrease  at  boundary.  If  the  boundary  pro¬ 
perties  are  not  different  from  those  in  the  bulk,  i.e.  n  =  n  . 

and  we  neglect  the  effect  of  the  self  magnetic  field,  the 

steady  state  solution  is  n  =  n  .  Thus  we  have  a  constant 

s 

density.  All  previous  authors  have  applied  only  the  first 
boundary  condition,  in  their  result,  the  coefficient  of  the 
Bessel  function  of  density  is  arbitrary.  By  the  second  bound¬ 
ary  condition  we  imposed,  the  coefficient  of  the  Bessel 
function  of  density  is  not  an  arbitrary  constant  but  depends 
on  the  density  at  the  walls .  We  believe  that  the  second  con¬ 
dition  is  necessary,  therefore  n  ^  n. 

s 

Equations  (48),  (51)  a  (5*0  and  (55)  are  derived  on 
the  assumption  that  the  use  of  the  relations  (46)  and  (47) 
and  the  assumption  y^Boz  <<  1  is  valid.  The  general  form  of 
equations  (48)  and  (51)  is  the  same  for  both  the  weak  field 
case  and  the  strong  field  case.  The  only  difference  is  in 
the  E0  dependence  of  y  and  T+  on  the  left  hand  side  of  (48). 

For  the  frequency  of  rotation  of  the  helical  wave,  equations  (54) 
and  (55)  are  obtained  for  weak  and  strong  field  cases 
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respectively.  It  is  possible  to  compute  the  critical  values 

of  E o  ,  k  and  a> i  from  equations  (48),  (51),  (54)  and  (55) 

for  the  two  extreme  field  cases . 

For  m  =  Oj  in  (48),  the  magnetic  field  does  not 

enter  the  expression  at  all.  A  solution  of  the  form  (28)  with 

m  =  0  is  a  pure  density  wave  with  no  screw  motion.  It  can 

easily  be  proved  from  equation  (48)  that  m  =  0  mode  will 

always  decay.  Since  any  perturbations  will  be  opposed  by 

diffusion,  neglecting  the  $  terms  corresponds  to  assuming 

that  the  force  due  to  the  induced  magnetic  field  of  the 

current  is  smaller  than  the  diffusion  force  which  maintains 

the  steady  state.  If  we  take  $  terms  into  account,  there 

will  be  some  values  of  $  ^  $  ,  the  critical  induced  magnetic 

c 

field,  and  it  causes  a  driving  force  larger  than  the  diffusion 

force,  instabilities  may  occur. 

The  helical  type  instability  which  corresponds  to 

m  =  ±1  are  considered  in  the  following.  The  threshold  for  the 

onset  of  instability  will  occur  when  (48)  becomes  an  equility, 

i.e.  o)2  =  0.  The  parameters  of  the  system  are  of  course  the 

wave  number  and  the  electric  and  magnetic  fields  applied  to 

the  plasma.  For  any  fixed  value  Z?(k)  obtained  from  (51), 

there  will  be  corresponding  values  of  E0  and  B0  which 

z  z 

will  represent  the  onset  of  growing  of  that  wave  number.  The 
general  behaviour  is  as  expected,  namely,  larger  magnetic 
fields  are  required  to  cause  instabilities  when  lower  electric 
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fields  are  applied. 

For  small  values  of  B0  ,  Zi(k)  approaches  a  constant 

value,  and  frequency  varies  linearly  with  B0  as  well  as 

z 

linearly  in  $ .  If  y_  =  y  ,  the  frequency  of  the  helix  yield  a 

minimum  frequency  of  zero,  there  will  be  no  helical  mode.  It 

is  to  be  expected  from  symmetry  that  either  direction  of 

propagation  of  the  helix  is  equally  possible.  For  y  <  y  , 

the  frequency  changes  sign,  which  indicates  a  change  in 

direction  of  the  propagating  waves. 

The  general  form  of  equations  (48),  (51),  (54)  and 

(55)  are  similar  to  the  corresponding  equations  of  other 

authors9510,  B0  «  - v -  rr^—  and  uii  «  (y  -  y,)Bo  ,  except 

z  y_  +  y+  Hjoz  -  +  z 

that  our  equations  contain  the  influence  of  azimuthal  magnetic 
field. 

In  order  to  get  approximate  solutions  of  the  above 

quantities  without  the  use  of  computers,  we  make  the  approxi- 
s 

mation  that  =-=r  <<  1  and  expand  the  Bessel  functions  in  terms 

bt, 

of  their  arguments  and  keep  terms  up  to  third  power  of  the 
arguments.  By  using  the  values  s  =  30cm/sec.,  =  l400y  sec 
and  R  =  0.058cm.  from  Hurwitz  and  McWhorter9,  the  approximation 
—  <<  1  seems  reasonable.  Using  the  expansions  of  the  Bessel 

Kt, 

function,  approximate  values  of  all  the  H’s  in  (48),  (51),  (54) 
and  (55)  are  evaluated.  Threshold  values  of  E0  B0  ,  wave 

Z  ! Z 

number  k  and  the  frequency  of  the  helix  may  then  be  obtained 
from  the  corresponding  equation. 
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3* 

Approximate  solutions  for  b(E0  )B0  and  Zf(k)  at 

z  z 

the  threshold  take  the  same  form  for  both  weak  and  strong 
field  cases.  The  threshold  value  for  Z?(k)  is  obtained  from 
equation  (51),  which  gives 


Z^Oc)  =  —■  (  Y(jKo R2  -  |6»r2r|  -T-2-  ^ 


^  1  +  §KgR2 

m  24^  ,  n-n 

3oR  Rf  -i— 
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64 


3x ( 24 ) 
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p  2  p  2  S  S 

°  RC  ~ 


■( TljK ° R  2  - 


9R2r2_s  n  nS'  1  1_R 
te°R  R£  -  ^  J  B0„* 


n 


(57) 

is 

Using  the  value  of  Zf(k)  in  (57),  we  get  also  an  approximate 

threshold  value  of  Eo  Bo 

z  z 


b(y_+y+)Boz= 


64x(2) 


4m6oR2pT  ^r-5-J,0x(i,8),s(T|K§R2-  |g§R2p|  -p-5-)3*  *  R 


n 


n 


Bo 
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Z^(k)  increases  as  the  surface  concentration  increases,  and 

the  threshold  for  instability  increases  to  larger  values  of 

E0  and  B  0 „ . 
z  z 

The  critical  frequency  a>i  is  obtained  from  (54)  for 

weak  electric  field  case 

1  Cn 

bi  =  2“(y_  “  y+)  q—  >  (59) 

d 

and  from  (55)  for  strong  electric  field  case 


y_y ,  C 

b  i  =  ( y  -  y,)(l  +  - ) 

yf.  +  y*  Cd 


(60) 


where 
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c  = 
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m  |-3oR2 
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R£  - 
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18  ,  3  t/-  2  p  2 
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n  n-n 

93oR2  si  s 
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^  n 


)  R  $ 


(61) 

Equation  (57)  gives  just  the  magnitude  of  k.  For 
a  wave  of  the  form  exp[i(a)t  +  m0  +  kz)],  it  is  evident  that 
only  the  relative  signs  of  k,  m  and  a)  are  of  significant. 

We  choose  the  signs  of  k  and  m  to  make  oi  >  0  in  (  59)  and  (  60) 
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and  to  satisfy  (58)  .  The  resulting  signs  of  k,  m  and  go  for 

various  directions  of  Eo  and  B0  are  the  same  as  those 

z  z 

given  in  Table  1.  of  Hurwitz  and  McWhorter.  The  sense  and 
motion  of  the  helical  wave  may  be  deduced  from  the  signs  of 
k ,  m  and  co . 
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